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We investigate the consequence of the dimension reduction on the magnetic anisotropy of FePt 
and CoPt nanoparticles. Using an extension of the magnetic anisotropy model of Neel, we show 
that, due to a statistical finite size effect, chemically disordered clusters can display a magnetic 
anisotropy energy (MAE) as high as 0.5 x lO'' J/m^, more than one order of magnitude higher than 
the bulk MAE. Concerning Llo ordered clusters, we show that the surface induces a reduction of the 
MAE as compared to the bulk, due to the symmetry breaking at the cluster surface, which modifies 
the chemical order. 
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I. INTRODUCTION 

The understanding and the control of the nanostruc- 
ture magnetic properties is a challenging field of re- 
search, for both fundamental and applied physics. At 
the nanoscale, finite size effects and surface contribu- 
tions become dominant and lead to a modification of 
the properties as compared to the bulk, such as the en- 
hancement of the orbital magnetic moment and of the 
magnetic anisotropy energy (MAE) at the less coordi- 
nated atoms. ^""^ However, technological applications of 
such structures for ultrahigh density magnetic storage 
media are limited by thermal fluctuations of the magne- 
tization direction. To overcome this so-called superpara- 
magnetic limit, intense efforts are devoted to the inves- 
tigation of materials with a high MAE. In such a view 
FePt^"— and CoPliii"— have become some of the most 
studied materials as they display a very large magne- 
tocrystalline anisotropy (over 5 x 10^ J/m^) when ordered 
in the LIq phase. 

The effect of the Llo ordering in FePt and CoPt on 
the magnetic anisotropy in bulk crystals or thin films 
has been widely studied on both experimental and the- 
oretical points of viewji^"— In the case of nanostruc- 
tures, as the surface modifies the local chemical environ- 
ment, it is expected to have potentially a large impact 
on the MAE. Although many experimental studies had 
been performed^"—, it is only recently that a clear im- 
pact of the dimension reduction on the MAE has been 
experimentally evidenced for CoPt nanoparticles.^^ By 
using model samples, and combining various techniques, 
the authors have indeed been able to reliably determine 
the MAE of CoPt clusters with a diameter around 3 nm: 
while the anisotropy of chemically disordered particles 
was found to be higher than that of the bulk phase, it 
was the opposite for chemically ordered particles. 

In this paper, we use an extension of the phenomeno- 
logical anisotropy model of Neel to investigate the effect 
of chemical ordering and size reduction in both FePt and 



CoPt clusters. We first show that this simple model is in 
good agreement with bulk and thin film results. Then, 
taking advantage of its simplicity, which allows us to per- 
form an extensive study of FePt and CoPt nanoparticles, 
we use this model to investigate the MAE of clusters. 
Taking into account finite size effects by mean of a sta- 
tistical analysis on a large number of nanostructures, we 
demonstrate the impact of LIq chemical ordering on the 
MAE and show how the surface and the finite size mod- 
ify the results as compared to the bulk. Contrary to the 
results on pure Fe and Co clusters^i^^, here the surface 
induces a lowering of the MAE for Llg ordered clusters. 

Using a simple phenomenological pair-interaction 
model to calculate an effect as subtle as the magnetic 
anisotropy may seem a quite hazardous approach as com- 
pared to more elaborate models, like ab-initio calcula- 
tions. It is indeed well known that a theoretical determi- 
nation of the MAE is a hard task: because it corresponds 
to very slight variations in the total energy as a function 
of the magnetization direction, highly converged results 
are required. Moreover, there is nowadays no available 
derivation of the Neel model from first principles. Nev- 
ertheless, this model, which is based on symmetries, has 
been found to be in good agreement with experiments 
on magnetic nanostructures^SiM. In particular it can 
remarkably well account for the prevailing effect of low- 
coordinated atoms on the MAE. This let us think that 
this empirical model must be valid to a certain point, and 
that it can give reliable tendencies regarding the evolu- 
tion of the MAE in nanostructures with various param- 
eters (size, shape, configuration...). Besides, let us note 
that a Neel expression for the local anisotropy is also 
used in atomic scale magnetic simulations based on a 
Heisenberg-type model. Unfortunately, even it can 
be quite successfully compared to some ab-initio calcu- 
lations and experiments (see further in the paper), we 
cannot demonstrate the validity of the Neel model in the 
precise case under study: there is indeed no available ab- 
initio calculations of the MAE for comparable CoPt or 
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FePt nanoparticles. However, in order to study the effect 
of size reduction and of chemical ordering in particles of 
a few nanometers (which is precisely the size of inter- 
est for ultra- high density magnetic storage applications), 
there is no other method that can reasonably be used, 
and that would be more precise. We are then convinced 
that investigating such complex systems, which are out 
of reach for first-principle calculations for the moment, 
will give valuable physical insights. Knowing the results 
predicted by this simple and crude model will then en- 
able further comparison with more advanced calculations 
in the future. It is also important to check if the predic- 
tions of the Neel model are correct in the case of alloy 
nanoparticles, and we hope that our work will stimulate 
experimental studies. 

II. MAGNETIC ANISOTROPY OF A 
BIMETALLIC ALLOY 

The magnetocrystalline anisotropy (MCA) arises from 
the coupling of the magnetization with the crystal lattice 
through the spin-orbit interactioni^S. Then, for a given 
magnetic atom, it is directly linked to its local environ- 
ment. For a whole structure, it is the signature of the 
system symmetries. It can be described using a phe- 
nomenological pair model, first introduced by L. NeeP^. 
In its original form, for a single element crystal, the MCA 
energy is written as a pair interaction over the first neigh- 
bors. The magnetic anisotropy energy of a dimer is ex- 
pressed, using the lowest order term as2^ 

-Epair = -i (m ■ ef (1) 

where m is the magnetization unit vector and e the unit 
vector along the bond. The prefactor L is the anisotropy 
parameter, which has to be determined and is generally 
fitted on bulk experiments. 

To consider a given structure, we need two more hy- 
pothesis. First we suppose that the magnetization direc- 
tion is homogeneous over the whole structure (macrospin 
approximation), which is a quite suitable approximation 
for nanostructures with dimensions close to the exchange 
length)^ Second, we also suppose that the L parame- 
ter is the same in the bulk and at the cluster surface, 
and can be extracted from bulk measurements. Even 
if one could envisage a variation of the anisotropy pa- 
rameter with the different environment in a structure, in 
fact the crude model of Neel with a uniform L extracted 
from bulk measurements has already shown a good agree- 
ment (qualitative at least but also quantitative) for such 
systems i^i^i^i As the purpose of our study is the predic- 
tion of the main tendencies, we will make this assump- 
tion in order to restrict the number of parameters in the 
model. On the basis of these two hypothesis, the total 
MCA energy, i?a, is evaluated by performing a sum over 
all the nearest neighbor pairs. This leads to a quadratic 
form, that can be diagonalized into the orthonormal ba- 
sis set (ei, 62,63). The normalization condition on m 



further simplifies the expression to a biaxial form 

EJV = -Ki{m-eif + K2{ra-e2f + K (2) 

where V is the volume of the structure, Kx and K2 are 
positive constants and K a constant that can be omitted. 
Note that, in general cases, the new basis is different from 
the one defining the atomic lattice. The direction ei (re- 
spectively 62) corresponds to the easy (respectively hard) 
magnetization direction, with the anisotropy constant Ki 
(respectively ^^2)- This shows that this model, which is 
limited to second degree terms in m, always leads to a 
biaxial second order magnetic anisotropy displaying two 
degenerated energy minima. At the equilibrium (at zero 
temperature), the magnetization lies along ei and KiV 
corresponds to the minimum energy needed to switch the 
magnetization from one minimum to the other one. Note 
that the value of K2 does not influence this switching en- 
ergy. For a bulk cubic crystal, this model would give 
a zero anisotropy, as the second order expression is not 
compatible with the cubic symmetry (for a more accu- 
rate description of the MCA in this case, one should use 
higher order terms in eq. [1] see ref. HJ)- However, this 
crude model is suitable to describe surface anisotropy 
where the cubic symmetry is locally broken. 

To extend this model to an alloy containing two types 
of atom, named A and B in the following, three types 
of neighbor pair have to be considered!^ two between 
identical atoms and one between A and B atoms. We 
introduce three different anisotropy parameters: LaAi 
Lbb and Lab respectively for A- A, B-B and A-B pairs. 
The MCA energy is then evaluated by summing eq. [T] 
over all the nearest neighbors pairs, taking into account 
the chemical nature of the two atoms. 

This model allows us to describe the effect of LIq chem- 
ical ordering on the MCA. Neglecting any tetragonal dis- 
tortion, the atomic site lattice considered is face centered 
cubic (FCC), and we focus on equiatomic alloys. In this 
case, the LIq phase corresponds to an alternation of pure 
A and pure B planes in the [001] direction, leading to a 
chemical anisotropy along the c axis. Oppositely, in the 
chemically disordered phase named Al, the atoms are 
randomly distributed on the lattice sites. The magni- 
tude of the chemical order is usually quantified by the 
long-range order (LRO) parameter S, which varies be- 
tween (totally disordered alloy) and 1 (perfectly ordered 
alloy)— . Defining two sublattices labelled a and /3, cor- 
responding to successive planes in the [001] direction, S 
is given by 

S' = 2(na-l/2) (3) 

where is the fraction of a sites occupied by atoms A. 
Note that changing the sign of S is equivalent to invert- 
ing the convention for the choice of the a and /3 planes. 
As this convention is arbitrary, the sign is meaningless 
and only the absolute value of S has to be considered. 
Conversely, the probability p for each lattice site to be 
occupied by an atom of the correct type {Am a. a plane 
and i? in a /3 plane) is directly related to 5 as p = . 
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To calculate the MCA of a bulk crystal of given or- 
der parameter, we consider the probability for a nearest 
neighbors pair to be a A-A, B-B or A-B bond and we 
calculate its average magnetic anisotropy parameter (i). 
Three different nearest neighbors pairs with equivalent 
(i) are obtained; those lying in the a planes, those ly- 
ing in the /3 planes and those linking atoms from two 
different planes. The corresponding {L) parameters are 
respectively Lq, Lp, and Lap: 

La = p^Laa + {1-pTLbb + '^p{^-p)Lab 
Lp = [l - p)'^Laa+P^Lbb + '2'P0- - p)Lab 
Lap = p{1-p)Laa+p{1-p)Lbb + [p^ + {'^-pf']LAB■ 
Yo^: a chemically disordered crystal {p = 1/2), the 
probability to have a A or _B atom on a given site is 
identical, and La — Lp ~ Lap'- the cubic symmetry is 
recovered and the MCA is zero. As soon as a chemical 
order arises in the crystal {S > 0), the cubic symmetry 
is broken and the second order MCA term has a finite 
value. Summing the different bond contributions in a 
FCC unit cell of volume Vccw, the MAE is obtained as 

^^a,cell — 2 {La + Lp — 2Lap) (m • Bc)^ 

= 2{LAA + LBB-2LAB)S\xn-e,f, (4) 

with Gc the unit vector aligned along the c axis of the 
crystal. It corresponds to a uniaxial anisotropy, along the 
[001] direction (c axis), and coincides with the chemical 
anisotropy direction. The anisotropy constant for the 
bulk alloy is then given by 

K{S) = -2{Laa + Lbb - 2L^b)5VKcU (5) 

and has a simple quadratic form. This evolution can be 
understood looking back at our model: as a pair inter- 
action is considered, the model involves products of the 
occupation probabilities for a given site, which are linear 
in S. If 2Lab > Laa + Lbb, the anisotropy constant is 
positive, meaning that the c axis is the easy magnetiza- 
tion axis with the anisotropy constant Ki — K{S). 

We now apply this model to the case of XPt crystals, 
where X = Fe or Co. In these alloys, a strong second 
order MCA accompanies the chemical ordering into the 
LIq phase, and originates both from a local chemical en- 
vironment change and from a tetragonal distortion of 
the FCC lattice. However, this slight deformation has 
a much weaker influence on the MCA^^ and we neglect 
it in the following, as it was supposed in the previous 
development. In order to apply our model, we need to 
determine the parameters Lxx, ^PtPt and Lxpt- As the 
magnetic moment on the Pt atoms is small in these crys- 
tals^' we neglect the magnetic anisotropy of the Pt-Pt 
bonds by taking Lptpt = 0. A straightforward relation 
is then obtained between Lxx and Lxpt by fitting equa- 
tion[5]on the experimental MAE value of the LIq ordered 
crystal. Consequently, there only remains one adjustable 
parameter and we finally estimate Lxx from measure- 
ments on pure samples using the original Neel model-- : 



taking into account the magnetostriction parameters, the 
model gives Lxx{r) = Lxx {ro) + {r-ro) (ro ) , where 
r and rg are respectively the bond lengths in the alloy 
and in the pure crystal. We have determined the param- 
eters for both FePt and CoPt (see table U. The Lxx 
parameter for X = Fc and X = Co has been calculated 
using the values given in ref.^"^ and a bond length of re- 
spectively 2.72 A for Fe and 2.71 A for Co. In order to 



XPt 


Klio 




Lxpt 




(J/m^) 


(meV) 


(meV) 


FePt 


7.10" 


0.23 


0.73 


CoPt 


5.10" 


0.91 


0.89 



TABLE I: Parameters Lxx and Lxpt of the Neel model for 
FePt and CoPt, determined as explained in the text from the 
Llo ordered bulk MAE and the bulk Fe and Co Neel model 
parameters-^. 



check the accuracy of our model, we compare some pre- 
dictions with bulk and thin film results. The variation of 
bulk FePt MAE versus S is compared to experiments^ 
and first principle calculations2i, and, as it can be seen 
from fig. [1] shows a good agreement. The model is also 
compared to thin film experimental results. The MAE of 
a Pt/Co/Pt(lll) multilayer with a Co thickness of one 
monolayer is 

Ea I f 3 \ 

= — SicoPt - T^icoco = 2Ks 6 

S Sat V 2 / 

where S is the film surface and is the surface atomic 
density. Ks is the usual surface anisotropy at the Co/Pt 
interface'^^ (the factor 2 on the right hand side of eq. [6] 
accounts for the two Co/Pt interfaces in the film). Us- 
ing the values of table U we find with our model Ks — 
1.17 mJ/m^, close to the value of 1.15 mJ/m^ reported 
in rel2^. 



III. CHEMICAL ORDER IN NANOPARTICLES 

In order to investigate the relationship between chem- 
ical order and magnetic anisotropy, we first need to pre- 
cisely describe chemical order in nanostructures. Indeed, 
due to the finite size (and consequently the surface), the 
description is not as simple as in the bulk. First, the 
numbers of a and /3 sites are not necessarily equal like in 
the bulk. The order parameter has to be adapted, so that 
a particle with a perfect alternation of A and B planes 
displays 5 = 1. Using fa the proportion of a sites, we 
use a more general definition of 5:^° 

5=tS^ (7) 

where before, the fraction of a sites occupied 

by a A atom, and x is the atomic concentration of A. 
Whereas for the bulk = 1/2 and eq. [3]is recovered, it 
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FIG. 1: Variation of the anisotropy constant K for a bulk 
FePt crystal versus the long-range chemical order parameter 
S. The results of our model (continuous line) are compared 
to experiments performed by S. Okamoto et al.^^ and to first- 
principle calculations performed by J.B. Staunton et al^ 

varies depending on the shape and the size of a nanopar- 
ticle. Moreover, note that inverting the convention for a 
and /? planes changes fa (for example, the truncated oc- 
tahedron cluster, as the one represented in fig. [2j would 
display = 0.507 or 0.493 depending on the conven- 
tion), but of course do not change the absolute value of 
S. Besides, a physical inversion of the atomic type on 
each lattice site will change S into —S, meaning that 
the physical chemical order is unchanged. However, in 
this case, the inversion leads to a cluster with a differ- 
ent stacking (in particular, the surface atoms are not the 
same) and possibly with different magnetic properties, as 
it will be discussed later. 

To obtain a complementary description of the chemical 
order we can also use a short-range order (SRO) parame- 
ter, since for small size systems short-range ordering may 
dominate and control the nanoparticle properties. Al- 
though the SRO and LRO descriptions are equivalent for 
an infinite bulk (where there exists a direct relation be- 
tween p, the probability for each site to be occupied by an 
atom of the correct type, SRO and LRO parameters^) it 
is not true in nanoparticles. As our anisotropy model is 
limited to the first neighbors, we also limit the SRO pa- 
rameter to the nearest neighbors and use the description 
of Cowley. '^^ With pi.m,n the probability for an atom i on 
a given site to have, in the (Z, m, n) direction, a nearest 
neighbor of the same kind, a local order parameter can 
be written as 

-J 1 Pl,m.n 

^Lm.n t 

X 

where x is, as before, the atomic concentration in A. In 
an infinite AB FCC crystal and for a given atom, 12 of 
these parameters can be calculated. However, regarding 
the Llo ordering, they are all equivalent (note that this 




FIG. 2: Sketch of a Llo ordered truncated octahedron clus- 
ter containing 2406 atoms. The two colors represent the two 
types of atoms. Note that the topmost and bottommost 
planes contain the same kind of atom. Consequently, two 
types of cluster are possible: either the dark atoms corre- 
spond to A-type atoms, or they correspond to B-type atoms. 
These two clusters display different fa factors, but have the 
same physical LRO parameter. 

is not true for all types of order): ai^rn,n — o-i foi' {I, m, n) 
directions in the (a, 6) plane and ai^rn,n = —o-i otherwise, 
with fli = 2p — 1 if atom i is A-type on a a site or B-type 
on a /? site and 1 — 2p otherwise. One can then 

define the SRO parameter a as the average of the over 
the whole structure. Taking into account the occupation 
probability of each site, it leads to: cr = {2p — 1)^. 

This parameter is directly linked to the LRO parame- 
ter a,s a — S^. Note that this quadratic relation is under- 
stood by the fact that the LRO description is based on 
the lattice sites occupation whereas the SRO description 
is based on nearest-neighbor pair correlations. 

In a nanoparticle, the SRO parameter has to be 
adapted due to the surface, as it was the case for S. 
For each atom i, we calculate a local order parameter ai 
as: 



where rii is the proportion of nearest-neighbors (NN) of 
the same chemical type as atom i, in the (a, b) plane, Xi it 
the proportion of NN of the same chemical type as atom 
i, among all its NN, and fi is the proportion of NN in the 
(a, b) plane among all the NN (if atom i is surrounded 
by 12 NN like in the bulk, fi = 1/3; if it is at a surface 
site, fi may vary). The SRO parameter a is then defined 
as the average of ai over the whole atoms. Of course, this 
definition of a is identical to the Cowley-type one for the 
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infinite bulk. 

An important feature, in the case of nanoparticles, is 
that the finite size allows us to explicitly take into ac- 
count the statistical distribution of configurations. In- 
deed, the atomic arrangement for a given ordering prob- 
ability p is not unique, resulting in a statistical distribu- 
tion of the chemical order parameters S and a. Then, 
the direct relation between S and p and the one between 
a and 5, established for an infinite system, are no more 
valid. 

In order to investigate the relation between p, a and 
S in nanoparticles, we have randomly generated a great 
number of nanoparticles with a given value of p. For each 
cluster, the true value of S and a is calculated and we 
then obtain a statistical distribution of chemical order 
parameters (and of cluster composition), related to the 
small atom number in the cluster. Note that the choice 
of the direction defining the a and /3 sublattices is com- 
pletely arbitrary, so that the [100] as well as the [010] 
or the [001] direction can be chosen to compute S and 
a. In order to have a physical description of the chemical 
anisotropy, we always choose the directions that gives, on 
one hand the highest value of j^l, and on the other hand, 
the highest value of a. For reasons explained further, 
we restrict our study to nanoparticles with a truncated 
octahedron shape. 

We first focus on the chemical order in clusters where 
each site has an equal probability to be occupied by a 
A or B atom (i.e. p = 1/2). For each cluster, we have 
calculated both LRO and SRO parameters. Contrary 
to the bulk where a unique S and a correspond to a 
given p, clusters display distributions of LRO and SRO, 
which are shown in fig. [3] As explained before, this is 
related to the small size of the particle and illustrates the 
distribution of configurations. This effect decreases when 
the diameter of the cluster increases, as the standard 
deviation of these parameters scales as the square root of 
the total number of atoms. Consequently, finite average 
values are found which means that the clusters cannot 
be considered as perfectly chemically disordered. The 
201 atoms clusters display (cr) = 0.037 and {\S\) = 0.094 
but note also that some clusters have a non-negligible 
probability to display chemical order parameter as high 
as CT = 0.08 or \S\ — 0.15. Surprisingly, there is no 
detectable correlation between LRO and SRO for these 
clusters. This means that, for these disordered clusters, 
these parameters probe the chemical order in different 
and complementary ways and none of them is found to 
be better than the other one. 

In partially ordered clusters {p ^ 1/2), the situation 
is quite different. Although a statistical distribution still 
exists, its width decreases with increasing p, which in- 
dicates that the number of configurations decreases. A 
clear correlation appear between LRO and SRO and 
their means follow the bulk relations: (S) « 2p — 1 and 
(a) « {Sf. _ 

In conclusion, we insist on the importance of the sta- 
tistical distribution of configuration for nanoparticles. 



(a) 



B 



- 201 atoms 

- 2406 atoms 



(b) 



0.00 



0.05 



0.10 



0.15 




— o — 201 atoms 

— » — 2406 atoms 



0.2 



0.3 



FIG. 3: Distribution of the SRO (a) and LRO (ISI) parame- 
ters for disordered clusters containing 201 atoms (circles) and 
2406 atoms (triangles). 



Clusters generated with a given p present a distribution 
of chemical order parameters, and conversely, a variety of 
configurations correspond to clusters with the same order 
parameter. This means that, for small particles, chemical 
order parameters are not sufficient to fully characterize 
the chemical arrangement of atoms in the cluster: it is 
impossible to sum up an entire chemical configuration 
with only one or two numerical parameters. A statisti- 
cal analysis is then essential to correctly investigate the 
link between chemical order and the physical properties 
in nanoparticles. 



IV. MAGNETIC ANISOTROPY OF 
NANOPARTICLES 

We now use this anisotropy model to investigate the 
consequence of the size reduction in FePt and CoPt 
nanoparticles, with respect to their chemical order. We 
focus on nearly spherical clusters, as it corresponds to 
the usual shape experimentally observed. Moreover, as 
the host lattice for the FePt and CoPt crystal is FCC, 
we have chosen to restrict our study to clusters with a 
perfect truncated octahedron (TO) shape (see figured]), 
which is the equilibrium shape in this case^i^ii^. The 
magnetic anisotropy of pure Co TO clusters has already 
been studied using the Neel model^^. It has been shown 
that, due to the symmetries of the TO shape, both sur- 
face and volume MCA are zero (the MCA of the facets 
compensates) and that a surface anisotropy only arises 
when additional facets are present. Therefore, the MAE 
of perfect TO clusters, studied with our model will ex- 
plicitly reveal the effects of size reduction and chemical 
ordering. Moreover, for the sake of simplicity, and in or- 
der to avoid any effect coming from the variation of the 
magnetic element concentration, we have only considered 
clusters with a composition around 50 % Pt. The dipo- 
lar interaction is also neglected in our study, since it is 
expected to play a significant role only for non-spherical 
clusters (shape anisotropy) 

Concerning the magnetic anisotropy, whereas the bulk 
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FePt or CoPt MCA is only uniaxial, the cluster finite size 
implies that this is no more true and the more general 
biaxial form (see eq. [2]) is always found. However, since 
the main point of interest is the minimum energy needed 
to switch the magnetization from an easy magnetization 
direction to the opposite one, we restrict our discussion 
to Ki. 

A. Chemically disordered nanoparticles 

We first focus on the case of fully disordered clusters 
{i.e. when p = 1/2, each site having an equal probability 
to be occupied by a X or a Ft atom). For each clus- 
ter, we have calculated Ki by summing eq. [T]over all the 
nearest neighbors bonds. We obtain a statistical distri- 
bution which is displayed on figure S] for different types of 
clusters. Just like the S and a distributions, the Ki dis- 
tribution is related to the small size of the cluster and its 
standard deviation scales as the square root of the total 
number of atoms in the cluster. Note also, that as for \S\ 
and cr, the mean value of Ki decreases when the cluster 
size increases, indicating that for a large number of atoms 
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FIG. 4: Distribution of the magnetic anisotropy Ki for FePt 
(close symbols) and CoPt (open symbols) disordered clusters 
containing 201 atoms - about 1.7 nm in diameter - (circles) 
and 2406 atoms - about 4.3 nm in diameter - (triangles). 

Besides, we have observed that the \S\ and Ki distri- 
butions are statistically independent which means that 
no detectable correlation exists between them. At first 
sight, this may seem surprising. In fact these two quan- 
tities are very different: S probes the long-range chem- 
ical order whereas Ki is sensitive to the local chemical 
order around the magnetic atoms. A perfectly LIq or- 
dered bulk structure with an antiphase boundary (i.e. 



an inversion in the X and Ft planes stacking sequence) 
exemplifies this difference: whereas S will be zero, the 
magnetic anisotropy will almost be unaffected. A small 
S is then not incompatible with a significant MCA. Note 
however that a strong correlation exists between \S\ and 
Ki when clusters display a higher chemical order, as it 
will be discussed later. 

The local character of the magnetic anisotropy let us 
expect a stronger link between Ki and the SRO param- 
eter, because they both depend on NN correlations. In- 
deed, we obtain a noticeable correlation between Ki and 
a. However, although the correlation coefficient between 
Ki and a is significantly higher than the one between 
Ki and S, and moreover increases with the particle size, 
it remains small. In particular, it corresponds to an in- 
significant variation of (Ki) with ct, as compared to the 
Ki statistical distribution (this can be seen in figure [7] 
for a ^ 0). In small disordered particles, the configu- 
ration distribution is responsible for a partial "decorre- 
lation" between the SRO and the magnetic anisotropy. 
Although the SRO parameter is more suited than the 
LRO one, it still fails to establish a definite link between 
chemical order and magnetic anisotropy. 

We have also investigated the spatial distribution of 
the cluster easy magnetization axis direction (ei vec- 
tor). Whereas the LIq ordering induces an anisotropy 
axis along the c axis, we do not find such a result for dis- 
ordered clusters. The ei direction distribution satisfies 
the cluster cubic symmetry and the three a, b and c axes 
are found to be equivalent. However, the distribution 
is inhomogeneous and a probability maximum occurs in 
the [110] direction (and the other equivalent directions), 
as shown in fig. [S) This means that the easy magnetiza- 
tion axis is rarely oriented towards the {001} or {111} 
facets of the TO cluster but generally points towards the 
edges connecting two {111} faces. Moreover, this direc- 
tion also corresponds to a maximum of the MAE: in the 
case of the 201 atom FeFt cluster presented in fig. [SJ the 
mean MAE is about 0.35 x 10^ J/m^ for an easy axis 
aligned along the [110] direction while it is about 0.20 
and 0.25 x 10^ J/m^ respectively in the [001] and [111] 
directions. 



B. Chemical ordering 

To study the chemical ordering effect in nanoparticles, 
clusters with different p values were generated, leading 
in each case to a finite width distribution for \S\ and 
a (except for p = I). In this way, a large number of 
configurations is obtained allowing us to investigate the 
relation between \S\, a and Ki. The results are repre- 
sented in fig. ini and [7] revealing the statistical nature of 
the link between S, a and Ki. These Figures are anal- 
ogous to the bulk K{S) curve but also shows the finite 
width distributions. 

Although in the bulk, simple relations have been found 
between Ki and S (quadratic variation) or a (linear vari- 
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FIG. 5: (color online) Statistical distribution of the easy mag- 
netization axis for a 201 atom FePt cluster, represented in 
a 3D view along the [111] direction. The color scale repre- 
sents the probability density for the easy magnetization axis 
to point in the corresponding direction, red (resp. blue) areas 
indicating the most (less) probable directions. 




|S| |S| 



FIG. 6: (color online) Variation of the anisotropy constant Ki 
versus \S\ for a 201 atoms FePt (a) and CoPt (b) cluster, and 
for a 2406 atoms FePt (c) and CoPt (d) cluster. The color 
scale represents in the (|5'|,/i'i) plane, the density of proba- 
bility for a cluster with chemical long range order parameter 
|S| to have a MCA Ki (the corresponding color bar is shown 
in inset). The dashed line corresponds to the bulk behavior. 



ation), the relation is more complex for clusters. For 
small \S\ or a, the mean Ki is higher than the one of the 
bulk and we find the result previously described: ver- 
sus S, the Ki distribution is constant (mean value and 
width) and versus tr, a very small increase of the mean 
Ki occurs. In this case, the magnetic anisotropy is dom- 
inated by the small size effect, in particular by the statis- 





FIG. 7: (color online) Variation of the anisotropy constant 
Ki versus a for a 201 atoms FePt (a) and CoPt (b) cluster. 
The color scale represents in the {a,Ki) plane, the density 
of probability for a cluster with chemical short range order 
parameter \a\ to have a MCA Ki (the corresponding color 
bar is shown in inset). The dashed line corresponds to the 
bulk behavior. 



tical distribution of the atoms on the lattice sites. Above 
a threshold {Sth or ath respectively for LRO and SRO 
parameters), (Ki) starts increasing and the distribution 
narrows (note that the threshold is the same for the two 
order parameters as att = S^/J. (Ki) is proportional to 
jS'P or to a, similarly to the bulk, which means that the 
chemical order is sufficient to induce a significant effect 
on the MAE. However, the increase is smaller than for 
the bulk and the perfectly ordered LIq cluster magnetic 
anisotropy appears to be weaker than for the bulk. As 
it will be discussed later, this is due to a surface effect. 
The case of CoPt clusters is even more complex, as the 
magnetic anisotropy displays two distinct limits for per- 
fectly ordered clusters. The difference is 0.32 x 10® J/m'^ 
for the 201 atoms CoPt clusters. This corresponds to two 
different LIq ordered clusters where the Co and Pt atoms 
have been exchanged. As previously discussed in section 
mil this does not changes the order parameter but, while 
in one case the top and bottom facets are made of Co 
atoms, in the other case they are made of Pt atoms, the 
magnetic properties are different (the MAE is higher in 
this case). This effect also exists for FePt clusters but is 
weaker (the energy difference is only 0.2 x 10* J/m'^) and 
it is not visible on the figures. Finally, as the cluster size 
increases, the threshold decreases and the Ki value for 
LIq ordered clusters increases, meaning that the behavior 
is getting closer to the bulk one. 

Our calculations have also revealed that, as the chemi- 
cal order increases, the easy magnetization direction gets 
closer to the c axis. When p is high enough (which cor- 
responds to \S\ above the threshold mentioned earlier), 
the easy axis is contained in a cone around the c axis 
whose summit angle decreases as p increases, as shown 
in figure [5] 

As a matter of fact, our results do not show any par- 
ticular difference between the description using a SRO or 
a LRO parameter. In fact, as long as the chemical or- 
der is significant, a coherence is developed over the whole 
particle (statistically, a significant LRO goes with a sig- 
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FIG. 8: Variation of the anisotropy cone angle ^ani (summit 
half angle) versus p for 201 and 2406 FePt atoms clusters. 
The anisotropy cone is defined as the one, around the c axis, 
containing 90% of the easy magnetization axes. 

nificant SRO) and both parameters are highly correlated 
with the magnetic anisotropy. The SRO description is 
found to be slightly more convenient only for disordered 
clusters. On the other hand, the theoretical variation of 
the magnetic anisotropy with the LRO parameter may 
be more useful for comparisons with experiments, or with 
more advanced calculations since the LRO pa- 
rameter is usually used to quantify the degree of order 
in CoPt and FePt alloys. Moreover, while there is to our 
knowledge no reported experimental measurement of a 
SRO parameter, for CoPt or FePt nanoparticles as those 
considered in our study, it is possible to determine their 
individual (using for instance transmission electron mi- 
croscopy) or their mean LRO parameter.—"— 



V. DISCUSSION 

To apprehend the joint effects of surface and of small 
size in XPt clusters (X — Co,Fe), we focus on two ex- 
treme cases: fully disordered and perfectly LIq ordered 
clusters. Theses systems show a completely different vari- 
ation of the MCA versus the cluster size, as it can be seen 
in figure [9] In the limit of large diameters, bulk results 
are recovered: disordered clusters have a zero MCA and 
Llo ordered clusters display MCA of 7 and 5 x 10^ J/m^, 
respectively for FePt and CoPt. For a finite diameter, 
chemically disordered (resp. LIq ordered) clusters dis- 
play a higher (resp. lower) MCA than the bulk one. 

For perfectly ordered clusters {S = 1), the variation of 
{Ki) versus the number N of atoms in a cluster can be 
adjusted by this expression 



(K,) [N] ^Kv + j^ (9) 

where Ky is the bulk MAE, at the corresponding 5, and 
the second term is a small size correction. This small 



size correction can be read in terms of a surface magnetic 
anisotropy because it scales with the surface to volume 
ratio T,/V . It could be alternatively written as i^TgE/y, 
with Ks the surface magnetic anisotropy. Accordingly, 
the total anisotropy energy for a particle can then be 
written, in an usual way, as the sum of a volume and a 
surface contribution: Ea — KyV + Ks^. However, in 
the present case Kg is negative, which is very unusual. 
The negative sign of Ks indicates a lowering of the MCA 
at the smallest sizes, as the chemical order is broken at 
the surface: due to missing neighbors, the contribution of 
the XPt bonds is reduced and the total MCA decreases. 
We find Ks = -0.42 and -0.33 mJ/m^ for the FePt 
clusters, and —0.32 and —0.23 mJ/m^ for the CoPt clus- 
ters. For each type of clusters, two values for Ka (and 
thus for Ks) are found due to the non-equivalence be- 
tween the two LIq stacking possibilities, which indicates 
that the surface effects depend on the type of facets: the 
clusters which have the highest anisotropy have pure Pt 
upper and lower [001] facets. This small size effect can- 
not be neglected for the smallest clusters: for less than 
2000 atoms (about 3.5 nm in diameter), the decrease is 
more than 10 % as compared to the bulk. It is worth to 
mention that the absolute value of this surface magnetic 
anisotropy is comparable to the highest one reported for 
pure clusters embedded in various matrices^"—. 

For chemically disordered clusters {S — 0) equation [5] 
cannot be used to account for the variation of (Ki) with 
N. Interestingly, in this case the evolution corresponds 
to the expression 



(Xi) {N) = Kv + j^ (10) 

indicating that in this case, the small size contribution 
does not originate from a surface effect. Here, the ef- 
fect scales with the inverse square root of the number of 
atoms, which is reminiscent of a statistical effect. Indeed, 
as discussed before, for small disordered clusters the sta- 
tistical distribution of the various chemical configurations 
is responsible for a concomitant (but uncorrelated) en- 
hancement of the mean long range order parameter and 
of the mean anisotropy constant. 

In the case of intermediate S values, the evolution of 
{Ki) can be adjusted by a compination of the two pre- 
ceeding expressions: 

(^i)W-^v- + ^ + ^, (11) 

where the two constants Ka and Ki, respectively account 
for the statistical/configuration effect and the surface ef- 
fect. The results of the fits are shown in table |lT] for 
different values of both for FePt and CoPt clusters. 
The relative weight of these two contributions gives us in- 
formation on the origin of the anisotropy variation with 
size. We find that the magnitude of the surface contri- 
bution (i.e. the absolute value of Kb) increases with the 
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chemical order parameter of the particles, while the be- 
havior of the statistical contribution is less intuitive: Ka 
is first increasing with S and then finally vanishes when 
S reaches 1. Note also that, since the surface and the 
configuration contributions to the MAE have opposites 
signs, there exists a range of S values where there is al- 
most no evolution of {Ki) with the particle size. 




lo' lo' lo' lo' 

Atom number 



FIG. 9: Variation of the mean anisotropy constant {K\) ver- 
sus the number of atoms in the cluster for FePt clusters of 
different l^l parameters. The dots correspond to the calcu- 
lated (A'l), the error bar representing the FWHM of the K\ 
distributions. The lines are obtained by fitting the calculated 
values with the expression proposed in the text. 
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0.0 


4.2 


0.0 


0.0 


3.0 


0.0 


0.5 


1.9 


6.4 


-5.9 


1.3 


4.0 


-3.8 


0.75 


4.0 


3.6 


-7.4 


2.8 


3.1 


-5.0 


1 


7.0 


0.0 


-9.2/-8.5 


5.0 


0.0 


-6.9/-5.1 



TABLE II: Volume anisotropy Kv and small size correction 
coefficients Ka and Kb, determined for several S using eg. llll 
from the results in fig [3 The two different values of Kb for 
Llo ordered clusters account for the two types of Llo stacking 
as discussed previously. Magnetic anisotropics are given in 
10^ J/m^ 

The predictions of the empirical Neel model can be 
of great help to analyze experimental magnetic mea- 
surements. Indeed, according to the present results, we 
can expect the chemically disordered particles to have a 
quite high anisotropy (much higher than that of the bulk 
phase), whereas it should be in fact impossible to obtain 
small chemically ordered nanoparticles having an intrin- 
sic MCA as high as that of the bulk LIq phase. This is 
consistent with our very recent experimental findings on 
CoPt nanoparticles In addition, our calculations show 
that the anisotropy constant is falling off rapidly when a 
particle is not fully ordered, which could also be responsi- 
ble of the experimental observation of a lower MAE than 



expected. This point may be particularly relevant, since 
it is predicted that the most favorable configurations for 
small particles correspond to S < 1 4^ We emphasize that 
the results reported in this paper have been obtained 
by completely letting aside any energetic consideration 
for the various chemical arrangements. Nevertheless, our 
model could easily be combined with a structural op- 
timization method, either to determine the MAE of a 
few particular structures predicted to be the most stable 
ones, or to determine the MAE distribution of a statis- 
tical ensemble of structures, taking into account their 
probability of existence. 

Another interesting feature predicted by our model is 
the significant MAE distribution, even for a given parti- 
cle size and a fixed chemical order parameter, which is 
due to the existence of a statistical set of chemical con- 
figurations. As discussed elsewhere,— this can be one of 
the major source of MAE distribution in realistic assem- 
blies of clusters (with a distribution of composition, size, 
shape...). This feature, which is a particularity of mag- 
netic alloys, should be directly visible in experimental 
measurements on well-defined samples. 



VI. CONCLUSION 

In this study, we have extended the anisotropy model 
of NeelSi to describe the magnetocrystalline anisotropy 
in mixed FePt and CoPt crystals. This phenomenological 
model allows us to describe the effect of Llo ordering on 
the MCA and appears to be in good agreement with pre- 
viously reported results on bulk and thin films. We have 
used this model to investigate the specific modifications 
occurring in nanosized clusters of FePt and CoPt with 
a truncated octahedron shape. In chemically disordered 
clusters, we have shown that there always exists a par- 
tial chemical order, due to the small number of atoms, 
and that, oppositely to the bulk, the MCA is non zero. 
On the contrary, for Lip ordered clusters, we have shown 
that surface induces a strong diminution of the magnetic 
anisotropy, due to the broken LIq periodicity at the clus- 
ter surface. 

Our finding points a lack in this type of ordered alloys 
as materials for high-density magnetic storage applica- 
tions, where a high switching energy KcsV is needed to 
stabilize the magnetization direction. In pure or chem- 
ically disordered clusters, the anisotropy constant in- 
creases when the size decreases, thanks to a positive sur- 
face MAE (see e.g. ref . [3113143^ . which attenuates the size 
variation of the switching energy. On the contrary, in 
the LIq ordered CoPt and FePt clusters, K^s decreases 
when the size decreases due to the negative surface MAE, 
which worthen the size variation of the switching energy. 
Here, surface is a drawback that may be compensated 
by playing on appropriate surface effects. Changing the 
surface anisotropy sign could be achieved using the high 
interface magnetic anisotropy that adds to surface mag- 
netic anisotropy when the cluster surface is in contact 
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with an appropriate non ferromagnetic material, like Pt^ 
or CoO"*®. Surface engineering, by embedding the clus- 
ters in the matrix or by synthesizing core-shell clusters 
with the non ferromagnetic shell around the magnetic 
core, appears to be complementary approaches to these 
LIq ordered materials, in order to overcome the super- 
paramagnetic limit in nanoscale clusters. 



Acknowledgments 



We gratefully acknowledge V. Repain and Y. Nahas 
(MPQ, Paris, France) for stimulating discussions. We 
acknowledge support from the European Community 
(STREP SFINx no. NMP2-CT-2003-505587). 



* Electronic address: "rohart@lps.u-psud.fr' 

^ P. Gambardella, S. Rusponi, M. Veronese, S. S. Dehsi, 
C. Grazioli, A. Dallmeyer, I. Cabria, R. Zeller, P. H. Ded- 
erichs, K. Kern, C. Carbone, and H. Brune, Science 300, 
1130 (2003). 

^ I. M. Billas, A. Chatelin, and W. A. de Heer, Science 265, 
1682 (1994). 

^ S. Rusponi, T. Cren, N. Weiss, M. Epple, P. Buluschek, 
L. Claude, and H. Brune, Nature Mater. 2, 546 (203). 

* M. Jamet, W. Wernsdorfer, C. Thirion, D. Mailly, 
V. Dupuis, P. Melinon, and A. Perez, Phys. Rev. Lett. 
86, 4676 (2001). 

^ S. Sun, C. B. Murray, D. Weller, L. Folks, and A. Moser, 

Science 287, 1989 (2000). 
^ Z. R. Dai, S. Sun, and Z. L. Wang, Nano Letters 1, 443 

(2001). 

B. Rellinghaus, S. Stappert, M. Acet, and E. F. Wasser- 
mann, J. Magn. Magn. Mat. 266, 142 (2003). 

® C. Antoniak, J. Lindner, M. Spasova, D. Sudfeld, M. Acet, 
M. Farle, K. Fauth, U. Wiedwald, H.-G. Boyen, P. Zie- 
mann, F. Wilhelm, A. Rogalev, and S. Sun, Phys. Rev. 
Lett. 97, 117201 (2006). 

^ C. Rong, D. Li, V. Nandwana, N. Poudyal, Y. Ding, 
Z. L. Wang, H. Zeng, and J. Ping Liu, Adv. Mater. 18, 
2984 (2006). 

^° C. Rong, N. Poudyal, G. S. Chaubey, V. Nandwana, 

R. Skomski, Y. Q. Wu, M. J. Kramer, and J. Ping Liu, 

J. Appl. Phys. 102, 43913 (2007). 
" A. C. C. Yu, M. Mizuno, Y. Sasaki, H. Kondo, and K. Hi- 

raga, Appl. Phys. Lett. 81, 3768 (2002). 
^2 Y. Xu, Z. G. Sun, Y. Qiang, and D. J. Sellmyer, J. Magn. 

Magn. Mat. 266, 164 (2003). 

C. Petit, S. Rusponi, and H. Brune, J. Appl. Phys. 95, 
4251 (2004). 

" X. Sun, 

Z. Y. Jia, Y. H. Huang, J. W. Harrell, D. E. Nikles, K. Sun, 

and L. M. Wang, J. Appl. Phys. 95, 6747 (2004). 

M. Mizuno, Y. Sasaki, M. Inoue, C. N. Chinnasamy, 

B. Jeyadevan, D. Hasegawa, T. Ogawa, M. Takahashi, 

K. Tohji, K. Sato, andS. Hisano, J. Appl. Phys. 97, 10J301 

(2005). 

L. Castaldi, K. Giannakopoulos, A. Travlos, D. Niarchos, 
S. Boukari, and E. Beaurepaire, J. Magn. Magn. Mat. 290- 
291, 544 (2005). 

F. Tournus, A. Tamion, N. Blanc, A. Hannour, L. Bar- 
dotti, B. Prevel, P. Ohresser, E. Bonet, T. Epicier, and 
V. Dupuis, Phys. Rev. B 77, 144411 (2008). 
W. Grange, I. Galanakis, M. Alouani, M. Maret, J. P. 
Kappler, and A. Rogalev, Phys. Rev. B 62, 1157 (2000). 
S. S. A. Razee, J. B. Staunton, B. Ginatempo, 
E. Bruno,and F.J. Pinski, Phys. Rev. B 64, 014411 (2001). 
^" O. Ersen, V. Parasote, V. Pierron-Bohnes, M. C. Cadeville, 



and C. Ulhaq-Bouillet, J. Appl. Phys. 93, 2987 (2003). 
J. B. Staunton, S. Ostanin, S. S. A. Razee, B. Gyorffy, 
L. Szunyogh, B. Ginatempo, and E. Bruno, J. Phys. Cond. 
Matt. 16, S5623 (2004). 

R. Skomski, J. Appl. Phys. 101, 9N517 (2007). 
M. Jamet, W. Wernsdorfer, C. Thirion, V. Dupuis, 
P. Melinon, A. Perez, and D. Mailly, Phys. Rev. B 69, 
024401 (2004). 

R. Morel, A. Brenac, C. Portemont, T. Deutsch, and 
L. Notin, J. Magn. Magn. Mat. 308, 296 (2007). 
D. A. Garanin, and H. Kachkachi, Phys. Rev. Lett. 90, 
065504 (2003); H. Kachkachi, and E. Bonet, Phys. Rev. B 
73, 224402 (2006). 
'^^ T. Jourdan, F. Langon, and A. Marty, Phys. Rev. B 75, 
094422 (2007). 

P. Bruno, Phys. Rev. B 39, 865 (1989). 
L. Neel, J. Phys. Radium 15, 225 (1954). 

A. Hubert and R. Schafer, Magnetic domains (Springer 
Berlin, 1998). 

3° F. G. Nix and W. Shockley Rev. Mod. Phys 10, 1 (1938) 
S. Okamoto, N. Kikuchi, O. Kitakami, T. Miyazaki, Y. Shi 
mada, and K. Fukamichi, Phys. Rev. B 66, 024413 (2002) 
N. McGee, M. Johnson, J. de Vries, and J. A. de Stegee 
J. Appl. Phys 73, 3418 (1993). 
J. M. Gowley, J. Appl. Phys 21, 24 (1950). 
R. van Hardeveld and F. Hartog, Surf. Sci. 15, 189 (1969) 
L. Favre, V. Dupuis, E. Bernstein, P. Melinon, A. Perez 
S. Stanescu, T. Epicier, J. -P. Simon, D. Babonneau, J.- 
M. Tonnerre, and J.-L. Hodeau, Phys. Rev. B 74, 014439 
(2006). 

^® We have checked this point by explicitly calculating 
the shape anisotropy. Indeed, because we are dealing 
with alloys, even the highly symmetric truncated octahe- 
dron shape does not give zero dipolar anisotropy (shape 
anisotropy), but according to our calculations it remains 
negligible. 

^'^ S. Stappert, B. ReUinghaus, M. Acet, and E. F. Wasser- 

mann, J. Cryst. Growth 252, 440 (2003). 

K. Sato, Y. Hirotsu, H. Mori, Z. Wang, and T. Hirayama, 

J. Appl. Phys. 98, 024308 (2005). 
^'^ R. V. Petrova, R. R. Vanfleet, D. R. Richardson, B. Yao, 

and K. R. Coffey, IEEE Trans. Magn. 41, 3202 (2005). 

S. Saita and S. Maenosono, J. Phys. Gond. Matt. 16, 6385 

(2004). 

M. Jamet, V. Dupuis, P. Melinon, G. Guiraud, A. Perez, 
W. Wernsdorfer, A. Traverse, and B. Baguenard, Phys. 
Rev. B 62, 493 (2000). 

M. Jamet, M. Negrier, V. Dupuis, J. Tuaillon-Combes, 
P. Melinon, A. Perez, W. Wernsdorfer, B. Barbara, and 

B. Baguenard, J. Mag. Mag. Mat. 237, 293 (2001). 

S. Rohart, G. Raufast, L. Favre, E. Bernstein, E. Bonet, 
and V. Dupuis, Phys. Rev. B 74, 104408 (2006). 



11 



M. MuUcr, P. Erhart, and K. Albe, Phys. Rev. B 76, 

155412 (2007). 

F. Tournus, S. Rohart, and V. Dupuis, IEEE Trans. Magn., 
in press. 



V. Skumryev, S. Stoyanov, Y. Zhang, G. Hadjipanayis, 
D. Givord, and J. Nogues, Nature 423, 850 (2003). 



